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Abstract 

The prediction and prevention of spinal injury is an important aspect of preventive health 
science. The spine, or vertebral column, represents a chain of 26 movable vertebral bodies, joint 
together by transversal viscoelastic intervertebral discs and longitudinal elastic tendons. This 
paper proposes a new locally-coupled loading-rate hypothesis, which states that the main cause of 
both soft- and hard-tissue spinal injury is a localized Euclidean jolt, or SE(3)— jolt, an impul- 
sive loading that strikes a localized spine in several coupled degrees-of-freedom simultaneously. 
To show this, based on the previously defined covariant force law, we formulate the coupled 
Newton-Euler dynamics of the local spinal motions and derive from it the corresponding coupled 
SE(3) — jolt dynamics. The SE(3) — jolt is the main cause of two basic forms of spinal injury: (i) 
hard-tissue injury of local translational dislocations; and (ii) soft-tissue injury of local rotational 
disclinations. Both the spinal dislocations and disclinations, as caused by the SE(3) — jolt, are 
described using the Cosserat multipolar viscoelastic continuum model. 

Keywords: localized spinal injury, coupled loading-rate hypothesis, coupled Newton-Euler 
dynamics, Euclidean jolt dynamics, spinal dislocations and disclinations 
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1 Introduction 



Normal function of the human spine is possible due to a complex interaction of its components (i.e., 
vertebrae, ligaments, discs, rib cage, and muscles). Age, trauma, spinal disorders, and a host of other 
parameters can disrupt this interaction to an extent that in certain cases surgery may be required 
to restore normal function. Several spinal disorders have been described in Go el et al 2006] from a 
mechanical perspective. An understanding of these disorders can assist in the design and development 
of spinal instrumentation. As biomechanics begins to be intertwined with tissue engineering, a better 
understanding of the particular disorders may also provide insight into 'biological' solutions. 

In particular, the center of rotation of the upper cervical spine is an important biomechanical 
landmark that is used to determine upper neck moment, particularly when evaluating injury risk in 
the automotive environment |Chancey et al 2007| . Also, new vehicle safety standards are designed 
to limit the amount of neck tension and extension seen by out-of-position motor vehicle occupants 
during airbag deployments. The criteria used to assess airbag injury risk are currently based on 
volunteer data and animal studies due to a lack of bending tolerance data for the adult cervical spine 
|Nightingale et al 2007| . 

Also, lumbar spine pathology accounts for billions of dollars in societal costs each year. Although 
the symptomatology of these conditions is relatively well understood, the mechanical changes in the 
spine are not. Previous direct measurements of lumbar spine mechanics have mostly been performed 
on cadavers. The methods for in vivo studies have included imaging, electrogoniometry, and mo- 
tion capture. Few studies have directly measured in vivo lumbar spine kinematics with in-dwelling 
bone pins. In vivo 3D motion of the entire lumbar spine has recently been tracked during gait in 
[Rozumalski ct al 2008 . Using a direct (pin-based) in vivo measurement method, the motion of the 
human lumbar spine during gait was found to be triaxial. This appears to be the first 3D motion anal- 
ysis of the entire lumbar spine using indwelling pins. The results were similar to previously published 
data derived from a variety of experimental methods. 

The traditional principal loading hypothesis [McElhaney and Myers 1993, Whiting and Zernicke 1998] , 
which describes general spinal injuries in terms of spinal tension, compression, bending, and shear, 
is insufficient to predict and prevent the cause of the back-pain syndrome. Its underlying mechanics 
is simply not accurate enough. On the other hand, to be recurrent, musculo-skeletal injury must 
be associated with a histological change, i.e., the modification of associated tissues within the body. 
However, incidences of functional musculoskeletal injury, e.g., lower back pain, generally shows little 
evidence of structural damage [Waddcll 1998 . The incidence of injury is likely to be a continuum 
ranging from little or no evidence of structural damage through to the observable damage of muscles, 
joints or bones. The changes underlying functional injuries are likely to consist of torn muscle fibers, 
stretched ligaments, subtle erosion of join tissues, and/or the application of pressure to nerves, all 
amounting to a disruption of function to varying degrees and a tendency toward spasm. 

For example, in a review of experimental studies on the role of mechanical stresses in the genesis 
of intervertebral disk degeneration and herniation [Rannou et al 2001] . the authors dismissed simple 
mechanical stimulations of functional vertebra as a cause of disk herniation, concluding instead that 
a complex mechanical stimulation combining forward and lateral bending of the spine followed by 
violent compression is needed to produce posterior herniation of the disk. Considering the use of 
models to estimate the risk of injury the authors emphasize the need to understand this complex 
interaction between the mechanical forces and the living body [Seidel and G riffin 2001. Compressive 
and shear loading increased significantly with exertion load, lifting velocity, and trunk asymmetry 
[Granata and Marras 1995] , Also, it has been stated that up to two-thirds of all back injuries have 
been associated with trunk rotation Kumar and Narayan 2006 . In addition, load-lifting in awkward 
environment places a person at risk for low back pain and injury [Reiser et al 2008] . These risks 
appear to be increased when facing up or down an inclined surface. 

The safe spinal motions (flexion/extension, lateral flexion and rotation) are governed by standard 
Euler's rotational intervertebral dynamics coupled to Newton's micro-translational dynamics. On the 
other hand, the unsafe spinal events, the main cause of spinal injuries, are caused by intervertebral 
SE(3)-jolts, the sharp and sudden, "delta"- (forces + torques) combined, localized both in time and 
in space. These localized intervertebral SE(3)-jolts do not belong to the standard Newton-Euler 
dynamics. The only way to monitor them would be to measure "in vivo" the rate of the combined 
(forces + torques)- rise. 
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It is well known that the mechanical properties of spinal ligaments and muscles are rate dependent. 
As elongation rate increases, ligaments generally exhibit higher stiffness, higher failure force, and 
smaller failure strain. Previous studies have shown that high-speed multiplanar loading causes soft 
tissue injury that is more severe as compared to sagittal loading. This paper proposes a new locally- 
coupled loading-rate hypothesis, which states that the main cause of both soft- and hard-tissue spinal 
injury is a localized Euclidean jolt, or SE(3) — jolt, an impulsive loading that strikes a localized spine 
in several coupled degrees-of-freedom (DOF) simultaneously. To show this, based on the previously 
defined covariant force law, we formulate the coupled Newton-Euler dynamics of the local spinal 
motions and derive from it the corresponding coupled SE(3) — jolt dynamics. The SE(3)— jolt is the 
main cause of two forms of local discontinuous spinal injury: (i) hard-tissue injury of local translational 
dislocations; and (ii) soft-tissue injury of local rotational disclinations. Both the spinal dislocations 
and disclinations, as caused by the SE(3) — jolt, are described using the Cosserat multipolar viscoelastic 
continuum model. 

While we can intuitively visualize the SE(3)-jolt, for the purpose of simulation we use the necessary 
simplified, decoupled approach (neglecting the 3D torque matrix and its coupling to the 3D force 
vector). Note that decoupling is a kind of linearization that prevents chaotic behavior, giving an 
illusion of full predictability. In this decoupled framework of reduced complexity, we define: 

The cause of hard spinal injuries (discus hernia) is a linear 3D-jolt vector hitting some interverte- 
bral joint - the time rate-of-change of a 3D-force vector (linear jolt = mass x linear jerk). 

The cause of soft spinal injuries (back-pain syndrome) is an angular 3-axial jolt hitting some 
intervertebral joint - the time rate-of-change of a 3-axial torque (angular jolt = inertia moment x 
angular jerk). 

This decoupled framework has been implemented in the Human Biodynamics Engine |I vance vic~ 2 005] , 
a world-class neuro-musculo-skeletal dynamics simulator (with 270 DOFs, the same number of 
equivalent muscular actuators and two-level neural reflex control), developed by the present au- 
thor at Defence Science and Technology Organization, Australia. This kinematically validated human 
motion simulator has been described in a series of papers and books [Ivancevic~an d Snos well 20011 
llvancevic" and Beagley 2003[ Ejin cevic 20021 1! vanccvic "20041 jlvancevic and Beagley 2005| , 
[Ivanccvic and Ivanccvic 2006a. Ivanccvic and Ivanccvic 2006b. Ivancevic and Ivancevic 2006c , 
Ivancevic and Ivancevic 2007c, Ivancevic 2006, Ivancevic and Ivancevic 2007a, Ivanccvic and Ivancevic 2006, 
llvancevic and Ivancevic 2007bl llvancevic and Ivancevic 2008] . 

2 The SE(3)— jolt: the main cause of spinal injury 

In the language of modern biodynamics [Ivancevic 2004, Ivancevic and Ivanccvic 2006a , 
[Ivancevic and Ivancevic 2006b[ llvancevic and Ivancevi c 2006c, llvancevic and Ivancevic 20 07d]. 
[Ivancevic and Ivancevic 2007ej . the general spinal motion is governed by the Euclidean SE(3)-group 
of 3D motions (see Figured]). Within the spinal SE(3)-group we have both SE(3)-kinematics (consist- 
ing of the spinal SE(3)-velocity and its two time derivatives: SE(3)-acceleration and SE(3)-jerk) and 
the spinal SE(3)-dynamics (consisting of SE(3)-momentum and its two time derivatives: SE(3)-force 
and SE(3)-jolt), which is the spinal kinematics x the spinal mass-inertia distribution. 

Informally, the locaUzed spinal SE(3)-jolE is a sharp and sudden change in the localized spinal 
SE(3)-force acting on the localized spinal mass-inertia distribution. That is, a 'delta'-change in a 3D 
force-vector coupled to a 3D torque-vector, striking the certain local point along the vertebral column. 
In other words, the localized spinal SE(3)-jolt is a sudden, sharp and discontinues shock in all 6 coupled 
dimensions of a local spinal point, within the three Cartesian (x, y, z)-translations and the three cor- 
responding Euler angles around the Cartesian axes: roll, pitch and yaw [Ivancevic and Bcagle y~2003| . 
If the SE(3)-jolt produces a mild shock to the spine, it causes mild, soft-tissue spinal injury, usually 

1 The mechanical SE(3)-jolt concept is based on the mathematical concept of higher— order tangency (rig- 
orously defined in terms of jet bundles of the head's configuration manifold) Ivancevic and Ivanccvic 2006c 
Ivanccvic and Ivancevic 2007c , as follows: When something hits the human head, or the head hits some external 
body, we have a collision. This is naturally described by the SE(3)-momentum, which is a nonlinear coupling of 3 linear 
Newtonian momenta with 3 angular Eulerian momenta. The tangent to the SE(3)— momentum, defined by the (absolute) 
time derivative, is the SE(3)— force. The second-order tangency is given by the SE(3)— jolt, which is the tangent to the 
SE(3)— force, also defined by the time derivative. 
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Figure 1: Human body representation in terms of SE(3)/SE(2)-groups of rigid-body motion, with 
the vertebral column represented as a chain of 26 flexibly-coupled SE(3)-groups. 



resulting in the back-pain sindrome. If the SE(3)-jolt produces a hard shock to the spine, it causes 
severe, hard-tissue spinal injury, with the total loss of movement. 

Therefore, we propose a new combined loading-rate hypothesis of the local spinal injury instead 
of the old principal loading hypothesis. This new hypothesis has actually been supported by a num- 
ber of individual studies, both experimental and numerical, as can be seen from the following brief 
review. One of the first dynamical studies of the head-neck system's response to impulsive loading 
was performed in |Misra and Chakravarty 1985 . The response of a human head/neck/torso system 
to shock was investigated in |Luo and Goldsmith 1991] . using a 3D numerical and physical models; 
the results indicated that the head, cervical muscles and disks in the lumbar region were subjected to 
the greatest force changes and thus were most likely to be injured. Time-dependent changes in the 
lumbar spine's resistance to bending was investigated in [Adams and Dolan 19 96] . with the objective 
to show how time-related factors might affect the risk of back injury; the results suggested that the 
risk of bending injury to the lumbar discs and ligaments would depend not only on the loads applied 
to the spine, but also on loading rate. Cyclic loading tests were performed by |Tsai et al 1998] to 
investigate the mechanical responses at different loading rates; the results indicated that faster load- 
ing rate generated greater stress decay, and disc herniation was more likely to occur under higher 
loading rate conditions. Anterior shear of spinal motion segments was experimentally investigated in 
|Yingling and McGill 1999 ; kinematics, kinetics, and resultant injuries were observed; dynamic load- 
ing and flexion of the specimens were found to increase the ultimate load at failure when compared 
with quasi-static loading and neutral postures. Experimental evidence concerning the distribution of 
forces and moments acting on the lumbar spine was reviewed in [Dola n~a,nd Ad ams 2001 , pointing 
out that it was necessary to distribute the overall forces and moments between (and within) different 
spinal structures, because it was the concentration of force which caused injury, and elicited pain. 
Small magnitudes of axial torque was shown to in [Drake et al 2005] to alter the failure mechanics 
of the intervertebral disc and vertebrae in combined loading situations. A finite element model of 
head and cervical spine based on the actual geometry of a human cadaver specimen was developed 
in [Zhang et al 2006 , which predicted the nonlinear moment-rotation relationship of human cervical 
spine. Vertebral end-plate fractures as a result of high-rate pressure loading were investigated in 
[Brown et al 200"8] , where a slightly exponential relationship was found between peak pressure and its 
rate of development. 

The localized spinal SE(3)-jolt is rigorously defined in terms of differential geometry 
[Ivancevic and Ivancevic 20 06c, Iva ncevic and Ivancevic 2007e] . Briefly, it is the absolute time-derivative 
of the covariant force 1-form (or, co-vector field) applied to the spine at a certain local point. With 
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this respect, recall that the fundamental law of biomechanics - the so-called covariant force law 
|Ivancevic and Ivancevic 2006b, I vancevic and Ivancevic 2006c[|Ivancevic and Ivancevic 200"7e] . states: 

Force co-vector field = Mass distribution x Acceleration vector-field, 

which is formally written (using the Einstein summation convention, with indices labelling the three 
local Cartesian translations and the corresponding three local Euler angles): 

F/j, = ra tiu a v ', (//, v = 1, 6 = 3 Cartesian + 3 Euler) 

where denotes the 6 covariant components of the localized spinal SE(3)-force co-vector field, 
Tn^u represents the 6x6 covariant components of the localized spinal inertia-metric tensor, while a v 
corresponds to the 6 contravariant components of localized spinal SE(3)-acceleration vector-field. 

Now, the covariant (absolute, Bianchi) time-derivative %{■) of the covariant SE(3)-force defines 
the corresponding localized spinal SE(3)-jolt co- vector field: 

j t {F,) = m^j t (a") = rTv, (6" + E>^ A ) , (1) 

where jjx(a>") denotes the 6 contravariant components of the localized spinal SE(3)-jerk vector-field 
and overdot (') denotes the time derivative. are the Christoffel's symbols of the Levi-Civita 

connection for the SE(3)-group, which are zero in case of pure Cartesian translations and nonzero in 
case of rotations as well as in the full-coupling of translations and rotations. 

In the following, we elaborate on the localized spinal SE(3)-jolt concept (using vector and ten- 
sor methods) and its biophysical consequences in the form of the localized spinal dislocations and 
disclinations. 



2.1 SE(3) — group of local spinal motions 

Briefly, the SE(3) — group of localized spinal motions is defined as a semidirect (noncommutative) 
product of 3D intervertebral rotations and 3D intervertebral micro-translations, 

SE(3) := SO(3) OR 3 . 

Its most important subgroups are the following (see Appendix for technical details): 



Subgroup 




Definition 




50(3), group of rotations 
in 3D (a spherical joint) 


Set of all proper orthogonal 
3 x 3 — rotational matrices 


SE(2), special Euclidean group 
in 2D (all planar motions) 


Se 


t of all 3 x 3 — matric 
cos 9 sin r x 
— sin 6 cos 6 r y 
1 


es: 


SO(2), group of rotations in 2D 
subgroup of 5-E(2)-group 
(a revolute joint) 


Set of all proper orthogonal 
2 x 2 — rotational matrices 
included in SE(2) — group 


M' 1 , group of translations in 3D 
(all spatial displacements) 


Euclidean 3D vector space 



In other words, the gauge SE(3) — group of intervertebral Euclidean micro-motions contains ma- 
R p 
1 



trices of the form 



where p is intervertebral 3D micro-translation vector and R is in- 



tervertebral 3D rotation matrix, given by the product R = R v ■ ■ Rg of the three Eulerian 
intervertebral rotations, roll = R„, pitch = R^, yaw = Rg, performed respectively about the 



x— axis by an angle if, about the y— axis by an angle ip, and about the z— axis by an angle 
[Ivancevic 20041 |Park and Chung 20051 llvancevic 2006] 1. 



(see 



Rip — 





COS If 

sin if 





— sin if 
cos If 



Rib — 



cos ip 


— sin ip 



sin %]) 

1 

cos ip 



Re = 



cos ( 
sin ( 




— sin 8 
cos 6 




Therefore, natural intervertebral SE(3)— dynamics is given by the coupling of Newtonian (trans- 
lational) and Eulerian (rotational) equations of intervertebral motion. 
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2.2 Localized spinal S E(3)— dynamics 

To support our locally-coupled loading-rate hypothesis, we formulate the coupled Newton-Euler 
dynamics of localized spinal motions within the SE(3)— group. The forced Newton-Euler equations 
read in vector (boldface) form 

Newton : p = Mv = F + p x w, (2) 
Euler : 7r=Icj = T + 7rXw + pXv, 

where x denotes the vector cross product!! 

M = Mij ; = diag{mi 1 m 2 ,m 3 } and I = I i3 ■= diag{h, I 2 , h}, (hj = 1,2,3) 

are spinal segment's (diagonal) mass and inertia matrices!! defining the localized spinal mass-inertia 
distribution, with principal inertia moments given in Cartesian coordinates (x, y, z) by volume integrals 

I\ = III p(z 2 + y 2 )dxdydz, I 2 = III p(x 2 + z 2 )dxdydz, I 3 = / / / p(x 2 + y 2 )dxdydz, 



dependent on localized spinal density p = p{x,y, 2), 

v = v l = [vi, v 2 , v 3 Y and id = uj 1 = [ui, uj 2 , w 3 ]* 

(where [ ]* denotes the vector transpose) are localized spinal linear and angular velocity vectorsQ (that 
is, column vectors), 

F = Fi — [F U F 2 , F 3 ] and T = T t = [T U T 2 , T 3 ] 

are gravitational and other external force and torque co-vectors (that is, row vectors) acting on the 
spine, 

p = p,. = Mv = \pi,p 2 ,P3] = [mivi,m 2 v 2 ,m 2 v 2 ] and 

7T EE TTi = IhJ = , 7T 2 , 7T 3 ] = [IlUJ\, I 2 U! 2 , I3OJ3] 

are localized spinal linear and angular momentum co-vectors. 
In tensor form, the forced Newton-Euler equations ([2]) read 

p % ee v j =Fi+ e J ikPj oj k , (i, j, k = 1, 2, 3) 
ni ee hjuji =T i + 4 k Tr j uj k + ei k p j v k , 

where the permutation symbol e\ k is defined as 

f+1 if k) is (1, 2, 3), (3, 1, 2) or (2, 3, 1), 
4 = | -1 if (*»j.*0 is (3, 2,1), (1,3, 2) or (2,1,3), 
I otherwise: i = j or j = k or k = i. 

In scalar form, the forced Newton-Euler equations expand as 

( p x = F\ — m 3 v 3 uj 2 + m 2 v 2 u 3 
Newton : < ^=^2 + "13^3^1 — miviuj 3 , (3) 



Euler 



Pi 


= F X 


— m 3 v 3 uj 2 + m 2 v 2 u> 3 




P-2 


= F 2 


+ m 3 v 3 uui — niiviuj 3 , 




P 3 


= F 3 


- m 2 v 2 ui + niiviuj 2 






= T X 


+ (m 2 - m 3 )v 2 v 3 + (I 2 


- I 3 V2W3 




= T 2 


+ (m 3 - mi)v-LV 3 + (I 3 


- h)uiuj 3 


K 


= T 3 


+ (mi - m 2 )v\v 2 + (h 


- I 2 )bJlL0 2 



2 Recall that the cross product u X v of two vectors u and v equals u X v = uvsiiiQxi, where 9 is the angle between 
u and v, while n is a unit vector perpendicular to the plane of u and v such that u and v form a right-handed system. 

3 In reality, mass and inertia matrices (M,I) are not diagonal but rather full 3x3 positive-definite symmetric 
matrices with coupled mass- and inertia-products. Even more realistic, fully-coupled mass— inertial properties of a 
spinal segment are defined by the single non-diagonal 6x6 positive— definite symmetric mass— inertia matrix Mse(3)i 
the so-called material metric tensor of the SE(3)— group, which has all nonzero mass— inertia coupling products. However, 
for simplicity, in this paper we shall consider only the simple case of two separate diagonal 3x3 matrices (M, I). 

4 In reality, u; is a 3 X 3 attitude matrix (see Appendix). However, for simplicity, we will stick to the (mostly) 
symmetrical translation— rotation vector form. 
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showing localized spinal mass and inertia couplings. 

Equations can be derived from the translational + rotational kinetic energy of the spine 

segment^ 

E k = iv'Mv + (4) 

or, in tensor form 

E = -M lj v i v j + -LjuftJ. 
2 3 2 3 

For this we use the Kirchhoff-Lagrangian equations (see, e.g., [Lamb 1932| ILeo nard 1997], or the 
original work of Kirchhoff in German) 

j t d v E k = 8 v E k x lu + F, (5) 

4[ d ^ E k = d u E k x lj + d v E k x v + T, 
at 

where d v E k = d^Ek — in tensor form these equations read 

j f d vi E = e j lk {d v] E)u k + F u 

j t d^E = el k (d^E)u k + e j lk (d v] E)v k +T l . 

Using (H])-©, localized spinal linear and angular momentum co- vectors are defined as 

p = d v E k , n = d u E k , 

or, in tensor form 

Pi = d v iE, TTi = d u iE, 
with their corresponding time derivatives, in vector form 

d ^ r d d a ^ 

P= dt P= dt dvE > ^dt^dt 9 ^ 

or, in tensor form 

d d a j? ■ d d 

Pi = dt P * = dt^ ni =dt ni = dt^ 

or, in scalar form 

P = [pi,P2,P3] = [miVi,m 2 V2,m 3 v 3 }, tt = [tti , 7r 2 , 7r 3 ] = [Ii£i>i, J 2 W2, h^]- 

While spinal healthy SE( 3) — dynamics is given by the coupled Newton-Euler micro-dynamics, the 
localized spinal injury is actually caused by the sharp and discontinuous change in this natural SE(3) 
micro-dynamics, in the form of the SE(3)— jolt, causing localized discontinuous spinal deformations, 
both translational dislocations and rotational disclinations. 

2.3 Localized spinal— injury dynamics: the SE(3)~ jolt 

The SE(3)— jolt, the actual cause of spinal injury (in the form of the localized spinal plastic deforma- 
tions), is defined as a coupled Newton+Euler jolt; in (co)vector form the SE(3)— jolt read^f] 



SE(3) - jolt 



Newton jolt : F = p p x u - pxtl>, 

Euler jolt : T = tt —it xw-?rxij-pxv-pxv, 



5 In a fully— coupled Newton— Euler localized spinal dynamics, instead of equation (J4j we would have spinal segment's 
kinetic energy defined by the inner product: 

Ek = - [pw \M S E(3) P w ] ■ 

6 Note that the derivative of the cross-product of two vectors follows the standard calculus product-rule: £ (u X v) = 
u X v + u X v. 
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where the linear and angular jolt co- vectors are 

F = Mv=[F u F 2 ,F 3 ], ± = 1w=[t 1 ,T2,t 3 ], 

where 

V = [i5i,U2,#3]*! w = [wi,i>2,w 3 ]*, 

are linear and angular jerk vectors. 

In tensor form, the SE(3)— jolt readfl 

*i - Pi ~ 4kPjU k - e{ kPj io k , (i, j, fe = 1, 2, 3) 

£ = -4k*j" k -4k^ k -4kPj yk -4ft-« fe ' 

in which the linear and angular jolt covectors are defined as 

F = Fi = Mv = My 5* = [Fi , F 2 , F 3 ] , 
f = T < =E5=Iyw J ' = [7i > r 2> T 3 ] J 

where v = ii 1 , and u — uj % are linear and angular jerk vectors. 
In scalar form, the SE(3) — jolt expands as 

pi - m 2 u}3V2 + m 3 (LU2V3 + V3OJ2) - 
p 2 + miuj^vi - m 3 u>ii 3 - m 3 v 3 ui + miUia; 3 , 
p3 - miW2«i + TO2W1W2 - v 2 oji - mividj 2 , 

7Ti - (m 2 - m 3 ) (« 3 U2 + - {h - h) (w 3 w 2 + ^2^3) , 

7T 2 + ("I! - m 3 ) (f 3 7ji + + (/1 - I 3 ) (CJ 3 W! + CJlW 3 ) , 

7T3 - (m-i - m-2) (w2«i + viv 2 ) - (h - I2) (w2^i + W1W2) • 

We remark here that the linear and angular momenta (p, 7r), forces (F, T) and jolts (F, T) are co- 
vectors (row vectors) , while the linear and angular velocities (v, uj) , accelerations (v, uj) and jerks (v, u>) 
are vectors (column vectors). This bio-physically means that the 'jerk' vector should not be confused 
with the 'jolt' co-vector. For example, the 'jerk' means shaking the head's own mass-inertia matrices 
(mainly in the atlanto-occipital and atlanto-axial joints), while the 'jolt'means actually hitting the 
head with some external mass-inertia matrices included in the 'hitting' SE(3)-jolt, or hitting some 
external static/massive body with the head (e.g., the ground - gravitational effect, or the wall - inertial 
effect). Consequently, the mass- less 'jerk' vector represents a (translational+rotational) non-collision 
effect that can cause only soft-tissue spinal injuries, while the inertial 'jolt' co-vector represents a 
(translational+rotational) collision effect that can cause hard-tissue spinal injuries. 

For example, while driving a car, the SE(3)-jerk of the head-neck system happens every time the 
driver brakes abruptly. On the other hand, the SE(3)-jolt means actual impact to the head. Similarly, 
the whiplash-jerk, caused by rear-end car collisions, is like a soft version of the high pitch-jolt caused 
by the boxing 'upper-cut'. Also, violently shaking the head left-right in the transverse plane is like a 
soft version of the high yaw-jolt caused by the boxing 'cross-cut'. 



Newton jolt 



Euler jolt 




2.4 Localized spinal dislocations and disclinations caused by the SE(3)— jolt 

Recall from introduction that for mild (soft-tissue) spinal injury, the best injury predictor is considered 
to be the product of localized spinal strain and strain rate, which is the standard isotropic viscoelastic 
continuum concept. To improve this standard concept, in this subsection, we consider spinal seg- 
ment (with a vertebral body, intervertebral disc and other visco-elastic tissue) as a 3D anisotropic 
multipolar Cosserat viscoelastic continuum [Cosserat an d Cosscrat 1898, Cosscrat and Cos scrat 1909, 
|Eringen 2002| , exhibiting coupled-stress-strain elastic properties. This non-standard continuum model 
is suitable for analyzing plastic (irreversible) deformations and fracture mechanics Bilby and Eshelby 1968 



7 In this paragraph the overdots actually denote the absolute Bianchi (covariant) time-derivative JTJ, so that the jolts 
retain the proper covector character, which would be lost if ordinary time derivatives are used. However, for the sake 
of simplicity and wider readability, we stick to the same overdot notation. 
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in multi-layered materials with microstructure (in which slips and bending of layers introduces ad- 
ditional degrees of freedom, non-existent in the standard continuum models; see [Mindlin 19651 
ILakes 1985] for physical characteristics and |Yang and Lakes 1981 Yang and Lakes 1 982|, 



[Park and Lakes 1986) for biomcchanical applications). 

The SE(3) — jolt (F, T) causes two types of localized spinal discontinuous deformations: 

1. The Newton jolt F can cause micro-translational dislocations, or discontinuities in the Cosserat 
translations; 

2. The Euler jolt T can cause micro-rotational disclinations, or discontinuities in the Cosserat 
rotations. 

For general treatment on dislocations and disclinations related to asymmetric discontinuous defor- 
mations in multipolar materials, see, e.g., |Jian and Xiao-ling 1 995, Yan g et al 200 1| . 

To precisely define localized spinal dislocations and disclinations, caused by the SE(3)— jolt (F, T), 
wc first define the coordinate co-frame, i.e., the set of basis 1 forms {rf.?''}, given in local coordi- 
nates x l — (x 1 , a; 2 , a: 3 ) = (x,y,z), attached to spinal segment's center-of-mass. Then, in the coor- 
dinate co-frame {dx 1 } we introduce the following set of spinal segment's plastic-deformation-related 
SE(3) — based differential p— forms (see [Ivancevic and Ivancevic 2006c[Hvancevic and Ivancevic 2007e] ): 

the dislocation current 1-form, J = Jidx 1 ; 

the dislocation density 2-form, a = \a.ij dx 1 A dx 3 ; 

the disclination current 2-form, S = \Sij dx 1 A dx 3 ; and 

the disclination density 3-form, Q = gj-Qyfc dx 1 A dx 3 A dx k , 
where A denotes the exterior wedge-product. According to Edelen Edelen 1980, Kadi c and Edelen 1 983J. 
these four SE(3)~ based differential forms satisfy the following set of continuity equations: 



a = dJ S, 


(6) 


Q = dS, 


(7) 


da = Q, 


(8) 


dQ = 0, 


(9) 



where d denotes the exterior derivative. 

In components, the simplest, fourth equation ([9]), representing the Bianchi identity, can be rewrit- 
ten as 

dQ = diQ[ijk] dx 1 A dx 1 A dx 3 A dx k = 0, 

where di = d/dx l , while denotes the skew-symmetric part of Oij.... 

Similarly, the third equation {5| in components reads 

■^Qijk dx 1 A dx 3 A dx k = SfcC*^] dx k A dx 1 A dx 3 , or 
Qijk = —6dkCX[ij] . 
The second equation ([7]) in components reads 

^Qijk dx 1 A dx 3 A dx k = -dkS^ dx k A dx 1 A dx 3 , or 

Qijk = 6dkS[ijy 

Finally, the first equation Q in components reads 



1 . , 1 

2 



6>ij dx 1 A dx 3 = {djJi — -Sij) dx 1 A dx 3 , or 



^QjJi Sij 



In words, we have: 



The 2-form equation ([B]) defines the time derivative a =hctij dx 1 Adx 3 of the dislocation density 
a as the (negative) sum of the disclination current S and the curl of the dislocation current J. 
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• The 3-form equation ([7]) states that the time derivative Q —-^Qijk dx % A dx J A dx k of the discli- 
nation density Q is the (negative) divergence of the disclination current S. 

• The 3-form equation ||5J) defines the disclination density Q as the divergence of the dislocation 
density a, that is, Q is the exact 3-form. 

• The Bianchi identity ([9]) follows from equation (J8j> by Poincare lemma [Ivancevic and Ivancevic 2 006c, 
llvancevic and Ivancevic 2007e] and states that the disclination density Q is conserved quantity, 
that is, Q is the closed 3-form. Also, every 4-form in 3D space is zero. 

From these equations, we can conclude that localized spinal dislocations and disclinations are 
mutually coupled by the underlaying SE(3) — group, which means that we cannot separately analyze 
translational and rotational spinal injuries — a fact which is not supported by the literature. 



3 Conclusion 

Based on the previously developed covariant force law, in this paper we have formulated a new 
coupled loading-rate hypothesis, which states that the main cause of localized spinal injury is an 
external SE(3) — jolt, an impulsive loading striking the spinal segment in several degrees-of-freedom, 
both rotational and translational, combined. To demonstrate this, we have developed the vector 
Newton-Euler mechanics on the Euclidean SE(3) — group of localized spinal micro-motions. In this 
way, we have precisely defined the concept of the SE(3) — jolt, which is a cause of rapid localized spinal 
discontinuous deformations: (i) mild rotational disclinations and (ii) severe translational dislocations. 
Based on the presented model, we argue that we cannot separately analyze localized spinal rotations 
from translations, as they are in reality coupled. To prevent spinal injuries we need to develop the 
internal SE(3)-jolt awareness. To maintain a healthy spine, we need to prevent localized SE(3)-jolts 
from striking any part of the spine in any human-motion or car-crash conditions. 



4 Appendix: The S E(3)— group 



Special Euclidean group SE(3) := SO (3) > R 3 , (the semidirect product of the group of rotations with 
the corresponding group of translations) , is the Lie group consisting of isometries of the Euclidean 3D 
space R 3 . 

An element of SE(3) is a pair (A, a) where A G 50(3) and a G R 3 . The action of SE{3) on R 3 is 
the rotation A followed by translation by the vector a and has the expression 

(A, a) ■ x = Ax + a. 

The Lie algebra of the Euclidean group SE(3) is se(3) = R 3 x R 3 with the Lie bracket 

[(£, u), (rj, v)] = (£ x rj, £ x v - r) x u). (10) 
Using homogeneous coordinates, we can represent SE(3) as follows, 



SE{3)= i ( n i ) S Gi(4,R) : R G 50(3), pei 3 



R p 
1 

with the action on R 3 given by the usual matrix-vector product when we identify R 3 with the section 
R 3 x {1} C R 4 . In particular, given 

R p 
1 



and q G R 3 , we have 

or as a matrix- vector product, 



G SE(3), 
■q = Rq+p, 



o i) (I)- rr 
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The Lie algebra of SE(3), denoted se(3), is given by 



se(3) 



U) v 





e M 4 (R) : w e so(3), t; £ 



where the attitude (or, angular velocity) matrix u> : M 3 — > so (3) is given by 






-w 2 









— U) : 










w = 

The exponential map, exp : se(3) — ► SE(3), is given by 

exp 

where 





V 







, 


)"( 



exp(w) Av 
1 



1 — cos ||oj|| ]HJ-sin|H 2 
A — I H -75 o> H 5 w , 



and exp(w) is given by the Rodriguez' formula, 



exp(w) = I + 



sin M 



1 — cos \\ui\ 



uj 2 . 
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